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Abstract 

The issue of the BRST symmetry in presence of the Gribov horizon is addressed in Euchdean 
Yang-Mills theories in the Landau gauge. The positivity of the Faddeev-Popov operator 
within the Gribov region enables us to convert the soft breaking of the BRST invariance ex- 
hibited by the Gribov- Zwanziger action into a non-local exact symmetry, displaying explicit 
dependence from the non-perturbative Gribov parameter. Despite its non-locality, this sym- 
metry turns out to be useful in order to establish non-perturbative Ward identities, allowing 
us to evaluate the vacuum expectation value of quantities which are BRST exact. These 
results are generalized to the refined Gribov- Zwanziger action introduced in |1J, which yields 
a gluon propagator which is non-vanishing at the origin in momentum space, and a ghost 
propagator which is not enhanced in the infrared. 
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1 Introduction 



The BRST symmetry is a fundamental tool for the study of gauge theories. This symmetry is 
at the origin of the Slavnov- Taylor identities which guarantee the renormalizability of Yang- Mills 
theories [2] . Furthermore, the use of the BRST charge enables us to identify the physical subspace 
of the theory, allowing to prove the unitarity of the S"— matrij|3. Though, reconciling this symme- 
try with the appearance of non-perturbative effects represents a non-trivial challenge in quantum 
field theory. 

In particular, the study of the BRST symmetry when the Gribov copies [3] are taken into account 
is of great relevance in order to unravel the non-perturbative features of confining Yang-Mills the- 
ories [U U] . While the Gribov copies can be neglected in the perturbative ultraviolet regime, they 
play a relevant role in the infrared, being closely related to gluon confinement. A partial resolution 
of the Gribov issue is provided by the Gribov-Zwanziger framework ^ El [7] , amounting to restrict 
the domain of integration in the functional integral to the Gribov region, whose boundary is the 
first Gribov horizon. Although the Gribov region is not completely free from Gribov copies |8], 
the restriction to this region can be achieved by adding to the Yang-Mills action a non-local term 
Sh, known as the horizon function [6l [7] 



Sh = 1^ J d^x h{x) , 

h{x) = g'r''A'^{M-'y'f'''Al, (1) 

where is the inverse of the Faddeev-Popov operator M"^ = -9^ (9^5"^ + gf^^'A"^) . The 

massive parameter 7 appearing in eq.([T]) is the Gribov parameter [3]. It is not a free parameter, 
being determined in a self-consistent way as a function of the gauge coupling constant g and of 
the invariant scale Aqcd through the so called horizon condition 

{h{x)) = A{N' - 1) . (2) 

Despite its non-locality, the horizon function can be cast in local form by means of a set of 
auxiliary fields, (jp^ , ,ZJlj' , tol^) , see Sect. 3. A local action is thus obtained from the Gribov 
horizon. Remarkably, it enjoys the important property of being renormalizable [6| [71 [9l [TOl [Tl ] [T2] . 

Having at our disposal a local renormalizable action encoding the restriction to the Gribov region, 
we can ask ourselves whether it possesses BRST invariance. Here, we need to be more precise 
about the question we are addressing. Till now, the BRST transformations we are referring to are 
those corresponding to the Faddeev-Popov Lagrangian with the inclusion of the auxiliary fields 
needed to localize the horizon function. As pointed out in [6l [7] , these fields give rise to a BRST 
quartet. Thus, we are considering the following BRST operator s: 



sAl = -Dfc^ , sc" = -gr'^'c^c' 



1 

2' 

sc"" = b\ 3^ = 0, 

„ — rCtb —rab „—:fit> 



— ab — ab — ab n to\ 



^Here, we are referring to theories for which the asymptotic fields can be consistently introduced. 
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where 6° is the Lagrange multipher enforcing the gauge condition and {if, c"") stand for the 
Faddeev-Popov ghosts. As a first answer to the previous question, one can check if the BRST 
transformations ([3]) leave the Gribov-Zwanziger action invariant. This is not the case. It turns out 
that the Gribov-Zwanziger action does not exhibit invariance under the BRST transformations 
([3]), which are broken by the presence of the horizoij^. Nevertheless, the resulting breaking term 
displays peculiar features. It is proportional to the Gribov parameter 7, being thus a soft break- 
ing. As such, it can be neglected in the deep ultraviolet region, where the notion of exact BRST 
symmetry is recovered. Moreover, it can be kept under control at the quantum level, due to the 
possibility of deriving generalized softly broken Slavnov- Taylor identities ensuring the renormal- 
izability of the Gribov-Zwanziger action. 

In [11], particular attention was devoted to the auxiliary fields {^1^, '^1^, cJ^*, cu"^) which, as much as 
the Faddeev-Popov ghosts (c", c"), develop their own dynamics, giving rise to non-perturbative ef- 
fects. We quote, for example, the dimension two condensate (O | i^^^{x)ip°i^{x) — LJ^{x)uj'jf{x)^ | O), 
easily written as the vacuum expectation value of a BRST exact quantity, i. e. 

(0 0) = (0 |. (x))| 0) . (4) 

Without the restriction to the Gribov region, such an expectation value would vanish. However, 
in the presence of the horizon, and thus for a non-vanishing Gribov parameter 7, it does not 
necessarily vanish. The explicit one-loop evaluation of expression (jl]) shows in fact that it is 
proportional to the Gribov parameter [TT] : 

(0 I ir^ixWti^) - I 0) = ^gVNiN' - 1)7^ . (5) 

That the vacuum expectation value of quantities which are BRST exact is non-vanishing has a 
simple understanding, due to the presence of a boundary in field space, namely, the Gribov hori- 
zon. We observe that the nilpotent BRST operator can be considered as an exterior derivative 
along the gauge orbit. However, even if the BRST exact quantity s {uj'lj'ip'^) can be seen as a 
total derivative, its integral, i.e. (O \ s {tJ^^f^) | O), picks up a boundary contribution, encoded in 
the explicit dependence of | {^'^'■p'^ — uj'I^uj'^'') | O) from the Gribov parameter. The existence 
of a soft breaking of the BRST symmetry plays an important role here. As noticed in \IT\, this 
breaking ensures that the parameter 7 is a physical parameter of the theory, entering thus the 
expression of the correlation functions of the theory. 

The analysis of the soft BRST breaking done in [11] has also revealed that the possibility of 
achieving an exact invariance through a suitable modification of the BRST operator by local 7- 
dependent terms has to be ruled out, due to the strong constraints imposed by dimensionality of 
the fields, Lorentz covariance and SU{N) color structure. Nevertheless, we cannot disregard the 
existence of a non-perturbative non-local symmetry of the Gribov-Zwanziger action, which would 
be compatible with the presence of the Gribov horizon. This is the issue we aim to face in this 
work. More specifically, we point out that the soft breaking of the BRST symmetry displayed by 
the Gribov-Zwanziger action can be converted into a non-local exact invariance. In other words, 
we shall be able to show that the BRST operator s and its soft breaking can be replaced by 
an operator s-y which corresponds to an exact invariance of the Gribov-Zwanziger action, while 

^ See Chapter V of [TT] for a detailed discussion on the soft breaking of the BRST symmetry. 
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exhibiting an explicit dependence from the Gribov parameter 7. Albeit the operator is non- 
local, it can be employed to investigate the dynamics of the auxiliary fields {Tp^^ , tp^^ ,uj^^ , uj^^ . 
We observe in fact that, imposing the condition 

(0|s^e(x)|0) = , (6) 

where 9(x) stands for a local field polynomial depending on the auxiliary fields v?^*, cJ^'', t^)^^) , 
will provide us a way to evaluate the BRST exact quantity (0 |sG(x)| 0). For example, from the 
condition (O (cJj^''(x)v9^''(x)) | O) = 0, we shall be able to obtain the expression of the condensate 
(0 I {^f{x)^f{x) - uf{x)uj''^\x)) I 0). Of course, the value obtained for (O | {^f^f - ujfujf) \ O) 
from condition ([6]) coincides with that which has been found by direct calculations [llj. 

We also underline that condition turns out to be compatible with the gap equation ([2]) de- 
termining the Gribov parameter 7. As we shall see, from the horizon condition ([2]) it follows 
that 

{0\s,{gr'^A;ix)ZJl%x))\0) = 0, (7) 
which is a particular case of condition ([6]). 

The paper is organized as follows. In Sect. 2, before facing the more complex case of the Gribov- 
Zwanziger action, we construct a toy model exhibiting soft breaking of the BRST invariance and 
reproducing features of the Gribov-Zwanziger action. This model will be employed to show how 
the operator can be introduced and how condition ([6]) enables us to evaluate the vacuum 
expectation value of BRST exact quantities. The Noether current corresponding to the exact 
s^-invariance will be derived. Sect. 3 is devoted to the Gribov-Zwanziger action. We introduce the 
operator and use it in order to evaluate the condensate (O K^^^'v^^'' — uJ^uj'jf') \ O). The rela- 
tionship between condition and the gap equation ([2]) is also discussed. In Sect. 4 we generalize 
the operator to the so called refined Gribov-Zwanziger model illj, in which the operator 
{^l^ip^ — cJ'lj'uj^^) is taken into account from the beginning. In Sect. 5 we collect our conclusion. 



2 A toy model and its BRST soft breaking 

Let us start by considering a scalar field (j) whose dynamics is described by the Euclidean action 

S, = l d'x Q0(-92)0 + ^0^) . (8) 

We look at a mechanism allowing us to modify in a smooth way the behavior of the correlation 
functions of the theory in the infrared, while leaving the ultraviolet behavior unmodified. As 
pointed out in [13] , this goal can be achieved by the mechanism of the soft breaking of the BRST 
symmetry. To that purpose, we introduce a BRST quartet (^, ip, tJ, oj) 

sip = u , su = , 

sU = ip , sTp = , (9) 

and 

s0 = , (10) 
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so that the operator s is nilpotent 

5^ = 0. (11) 

The fields {uJ,uj) are anti-commuting and have ghost number (—1,1), while are a pair of 

complex conjugate commuting fields, carrying vanishing ghost number. For the BRST invariant 
action we write 

>S'inv = Sff) + 5'oxact ; (12) 



where 



5'exact = S j d'^X [uj{-d'^)(f + p^UJif) 

d^x{i^{-d'^)^-uj{-d'^)uj + p^ij^'^-ujuj)) . (13) 



The action ^inv has the same physical content of S^f,. This follows by observing that the integration 
over the fields {Tp, (p, uJ, uj) amounts to introduce a unity factor in the path integral. 

Further, we introduce a soft breaking of the BRST invariance, through the term 

S^ = ^ [ d^x(l)iip-ip) , (14) 



2 „ 

where "(9 is a mass parameter, here introduced by hand, which will enable us to modify the large 
distance behavior of the correlation functions (0(xi)....0(x„)). Thus, for the starting action S, we 
have 

•S* = 5'inv + "S"^ , (15) 

and 

sS* = — / d^x(j)uj , (16) 

showing that the BRST symmetry is softly broken. It is interesting to note that the fields 
(j^,ip,Lj,u) can be integrated out also in the presence of the breaking term S^, yielding the 
nonlocal action 

The fields (jp, ip, ZJ, u) can be seen thus as auxiliary fields needed to localize the term ^_gf^^i^ 
in expression f|T7|) . In this sense, these fields play a role analogous to that of the auxiliary fields 



(JpI^ , ipl^ ,uJ^ , u'^^) needed to localize the horizon function of the Gribov-Zwanziger action. In the 
same way, the parameter ^9 is akin to the Gribov parameter 7. Despite the presence of the soft 
breaking, eq.f lTB]) . the action 5* turns out to be renormalizable. This is due to the fact that the 
BRST quartet {ip, (p, ZJ, u) is coupled in a linear way to the scalar field 0, as exhibited by expression 
(fT^ . As a consequence, the equations of motion of the fields (jp,ip,ZJ,uj) acquire the meaning of 
linearly broken Ward identities, namely 

^ = (-9 +P V-y0, 



6uj 



-(-9^+p> 
(-92 + p'p . 
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These identities imply that the most general counterterm turns out to be independent from the 
fields of the BRST quartet. It depends only on the scalar field 0, ensuring thus that the theory 
is renormalizable. In other words, the soft breaking of the BRST symmetry can be kept under 
control at the quantum level. 

Let us give a look at the propagators: 

mM-k)) 

From expressions (fT9l) it is apparent that the correlation functions of the scalar field, (0(a;i)....0(x„)), 
are modified in the infrared region, due to the presence of the soft parameters p). Moreover, 
one also observes that, in the absence of the BRST breaking, i.e. when = 0, the propagator of 
the scalar field turns out to be independent from the parameter p, which becomes an unphys- 
ical gauge parameter. As a consequence, in the absence of the breaking, the Green's functions 
{(f){xi)....(f){xn)) will be independent from p as well. However, for non-vanishing i}, the parameter 
p is no longer a gauge parameter, entering the expression of 



+ p^k^ + f ' 
1 

Yk^ + p^k^ + ^ ' 
1 

"TFTp2pTf ■ 

(19) 



The framework outlined here gives us a way to introduce infrared effects in a local and renor- 
malizable fashion. Nevertheless, the presence of the soft breaking of the BRST symmetry does 
not enable us to make use of the BRST operator s in order to characterize the vacuum of the 
theory. Due to the lack of a conserved charge, one cannot impose that the vacuum is annihilated 
by the BRST operator. This means that the vacuum expectation value of quantities which are 
BRST exact does not necessarily vanish, i.e. 

(0|se(x)|0) ^ . (20) 

So far, the only way to access correlation functions of the type of (!20l) is through explicit computa- 
tions. Let us give an example by considering the dimension two condensate {{'ip{x)ip{x) — ZJ{x)uj{x))) , 
which has precisely the form of (1201) . i.e. 

{{i^{x)<^{x) -uJ{x)uj{x))) = {s{uJ{x)<^{x))) . (21) 

This condensate can be seen as the analogue of the condensate (O | {^'jj'^l^ — TJ'^^ujI^) | O) occurring 
in the case of the Gribov-Zwanziger theory. Expression ( 12T1) can be obtained by differentiating 
the vacuum energy Ey with respect to the parameter p. In fact, from 

e"^^" = j [D<^] e-^ , (22) 

it follows that 

W = — / m] e-s — = ^^^^ - • ^22) 
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Making use of the dimensional regularization, the exphcit evaluation of the one-loop vacuum 
energy gives 

(24) 



{{Lfip - ujuj)) = {s (uip)) = -— / ,4 , , , o o . , ^ , , (25) 



^^ = \j^^ (log {p' + pV + y) - log ip' + P')) 

Thus, 

where the integral in the right hand side of eq. fl25|) is convergent. Therefore, we see that the con- 
densate {{^ip — ujuo)) is non- vanishing for non- vanishing -d. This example illustrates in a simple 
way that the we cannot access the vacuum state of the theory with the operator s. 

From expression (p5l) one also sees that the condensate ((^(/? — ujuo)) is non-vanishing when the 
parameter p is set to zero, i.e. 



r d^p 



iiW - ^^)),=o = -T / ,o 4 , o ^ . (26) 



Even if the operator {ifif — uJu) were not included in the starting action, it would have been 
generated dynamically, as it is apparent from eq.( l26i) . Therefore, the introduction of the parameter 
p enables us to take into account the operator {if(f — uJu) from the beginning. As we shall see, a 
similar situation will be encountered in the case of the Gribov-Zwanziger action. 



2.1 The modified BRST operator 

Let us look now at the possibility of finding a modified BRST operator which corresponds to an 
exact invariance of the action S. We shall rely on the property that the breaking term, eg. (1161) . is 
quadratic in the fields, due to the fact that the BRST quartet couples linearly to the scalar field 
0, see expression Let us consider the equation of motion of the field ZJ, i.e. 

-{-d^ + P^)lu{x) . (27) 



6uj{x 



We observe that the operator (—9^ + p^) is positive definite, so that its inverse does exist. As a 
consequence, eq.(l271) can be solved for uj{x), yielding 



where 



1 f d^p I 



[-d-'+P^U J (27r)>2 + p: 
Therefore, equation f|T6|) can be rewritten as 

f . 1 SS 
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showing that the BRST breaking of eg. (1161) can be cast in the form of a contact term, i.e. of a 
term related to the equations of motion. Equation (!30l) imphes that the action S is left invariant 
by the modified nilpotent operator 

§^ 1 



S^UJ 


= Tp + 




= 0, 




= ^ , 


S^UJ 


= 0, 




= 0, 




s^S = 



2 (-92 + p2) 



(31) 

so that 

. (32) 

The operator exhibits explicit dependence from the parameter Moreover, it reduces to the 
operator s when = 0. One should notice, however, that the operator is non-local. As such, it 
cannot be used to analyse the renormalizability properties of the model. Nevertheless, its existence 
is helpful in order to evaluate the expectation value of BRST exact quantities, by requiring that 

(0|s^e(x)|0) = . (33) 

To give an example of the usefulness of eq.(l33!l. let us compute again the condensate {(jpip — uJuj)), 
by making use of condition ( l33l) . From the equation 



it follows that 



{s^ {ZJ^)) = , (34) 



{{Tp^ - Uu)) + ^ ( ^-———^ ) = . (35) 



2 Yi-d^ + p^) 



Thus 



"Ty ^(27r)^(27r)>2 + p2 



^92 f dS 1 



2 J (27r)^g2 + p2 
Finally, from equations ( |T9l) we obtain 



W1M-1)} ■ (36) 



in agreement with expression (^^. We see thus that condition ( 1551) gives us a practical way of 
evaluating correlators of the kind (0 |s0(x)| 0). 
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2.2 The Noether current 



Let us conclude the analysis of the toy model by deriving the Noether current corresponding to 
the exact invariance (l32l) . To that purpose, we make use of the functional form of the operator 
s^, i.e. 

wix) = i^Tpix) + yA(x) 1 + ^(a^)^^ > (38) 

where A{x) stands for 



H^) = I Ar^j^-^'/'d/) = / Att^^— ^e-^^(^-^V(2/) , (39) 



so that 



i-dl + p')A{x) = 0(x) . (40) 

Acting with the operator w{x) on the action S one finds 

/ ^2 \ ^2 

w{x)S = i ifix) + y A(x) j {d'^u - p^uj) + {-d'^Tp + p'^Tp) oj + y 0u; 

= (Ipd^uj - (d^Tp) u;) + y {d^ (Ad^cu - (d^A) cu) + ( (5^ - p') A) cu) + y 0a; 

(41) 

Making use of eq- fjlUl) . it follows 

wix)S = d^3, , (42) 

where the Noether current j'^ is given by 

_ ^2 

if, = {(pdf,uj - {d^ip) ^^) + y {^d^,uj - (9^A) uj) . (43) 



2.3 Preliminary considerations 



Several remarks can be made from the example considered here. 



• The mechanism of the soft breaking of the BRST enables us to modify in a smooth way the 
correlation functions of the theory in the infrared, due to the presence of the soft parameters 
(^,P)- 

• The introduction of the BRST quartet and of the corresponding breaking term are done in 
a way which preserves locality and renormalizability. 
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• However, due to the presence of the breaking, the BRST operator s does not allow us to 
characterize the vacuum of the theory. Nevertheless, it turns out that the action S is left 
invariant by a modified operator s^, which exhibits explicit dependence from the parameter 
while reducing to the BRST operator s for vanishing ^9. Although non-local, the new 
operator can be employed to obtain the explicit expression of quantities which are BRST 
exact. 



3 The Gribov-Zwanziger action 

3.1 The non-local horizon function and its localization 

Let us start by giving a short overview of the Gribov-Zwanziger framework, enabling us to imple- 
ment the restriction in the Euclidean functional integral to the Gribov region Q [3] , defined as the 
set of field configurations fulfilling the Landau gauge condition and for which the Faddeev-Popov 
operator, 

M"' = -d^ {8^6'^' + gr'A^) , (44) 

is strictly positive, namely 

n = {A"^, d^A''^ = 0, M"^ > 0} . (45) 
This amounts to adding to the Yang-Mills action 



Syu = - I d'xF;,F^^^ , (46) 



the horizon term [6l [7] 

Sh = 1^ j d'^x h{x) , 
h{x) = g^r'^A'^iM-Yf^'A'',, (47) 
where the non-local kernel (Ai^^y^ stands for the inverse of the Faddeev-Popov operator 

M'^'ix) Um-\x, y))'"] = 6\x - y)6'' . (48) 
Thus, for the partition function implementing the restriction to Q one has [6l [7] 

/ DA 6{dA) detM e'^^'^' = [ DA 6{dA) detM e-(^^«+^'^) , 
Jq J 



(49) 



where the Gribov parameter 7 is determined by the horizon condition [6l [7] 

{h{x)) = A{N^ - 1) . (50) 

As already mentioned, the non-local horizon function can be localized through the introduction 
of a suitable set of additional fields [6l [7] 



e"^'^ = / DipDlpDiuDiJ e"""- , (51) 



j DipDlpDuDcJ e"^'° 
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with 



The fields {^1^, ^1^) are a pair of complex conjugate bosonic fields, while (p^, cu^^^ are anti- 

commuting fields. Each of these fields has 4 (A^^ ~ 1)^ components. Therefore, a local action is 
obtained from the restriction to the Gribov region, namely 

/ DA 6{dA) det M e~^^^^ = DA DipDipDuDuJ DbDcDc e~^«2 , (53) 
Jn J 

where Sqz is the Gribov-Zwanziger action, given by 

Sgz = Sym + j d^x {b^d^A"^ - c'^^W^^'c*) + j d^x TW^VjT + K'^^'^'^l') 

-j'g I d'x l^r'^A-,^'; + r^^A^if^ + ^ - ^) • 

The constant term, 4 (A^^ — 1) 7^, in expression (1541) is introduced in order to rewrite the horizon 
condition (ISUI) as 

11^ = 0. (55) 
with r being the effective action obtained from Sqz 



To the first order, the gap equation fl55|) takes the form 



Let us also observe that condition (1551) is equivalent to 

(0 0) + (0 \gr'''Alt,^ 0) = -8 [N' - l) f , (5^ 

which, owing to the discrete symmetry of the action ^gz 















- 





amn -rnc ^nc) ^ (59) 

becomes 

(0 \gr''Ay'^\ 0) = (0 l^^r'^A^^^ I 0) = -4 [N' - l) 7' = -7' Hx)) . (60) 

Equation fl60|) shows the relevance of the dimension two condensates \gf"'^'^A'^^Lp'^^\ O), \gf"'^'^A'^^'^^\ O), 
within the Gribov-Zwanziger approach. 
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3.2 Soft breaking of the BRST invariance 



As already underlined, the auxiliary fields {^'^, uJ^, tu^'') needed to localize the horizon func- 
tion (1471) give rise to a BRST quartet, meaning that they are assembled in two BRST doublets 
[HI [7], see eqs.(l3]). Acting thus with the BRST operator s on the action (15^ . one finds 

-^'g I d'x {c^d;^ ir"" {^'; + ^;)) + r'^Ay;) , (ei) 

which shows that the BRST symmetry is broken. One notes that the first term in the right-hand 
side of expression fl6T|) can be written as the BRST variation of gf^^"^ {^D^Jc^^ [djlflf) v^™^, i.e. 

gfatm ^f^g-ac^ ^n.c _ ^g-ac^ ^mc) ^ _^ ^^fabrn ^g-ac^ ^mc) ^ (g2) 

SO that it does not correspond to a real breaking. The terms which give rise to the true BRST soft 
breaking are those proportional to the Gribov parameter 7. The expression gf^^"^ [dI^c^) (^duZJl^) 
is usually introduced in the starting Gribov-Zwanziger action. It can be obtained through a change 
of variables [6i[7] in expression (!53|) . the corresponding Jacobian being field independent. However, 
for the purposes of the present work, we keep the term gf""^"^ (D^^c^) {^{djp^^^ (p^^ — {djj]°-^^ ^^'^) 
as it stands. 



3.3 The operator s 



To construct the generalized operator s^, we proceed as in the toy model considered in the previous 
section. From the Gribov-Zwanziger action, we get 



SSgz ^ ^abj,c 



(63) 



We observe now that, within the Gribov region fi, eq. fHSl) . the Faddeev- Popov operator is 
strictly positive, so that its inverse (A^~^)" does exist. This property enables us to solve equations 
(1631) for the fields and uJ'^^ obtaining 



c\x\ 



be 



5S, 



GZ 



i-^da SSgz 



dS \{M-\x,y)) 



da 



SSgz 



(64) 



Thus, equation fl6T]) can be rewritten as 

sSgz = j d'x (^-{D^K). [M'^) 



lymc SSgz 



bm SSgz 

xy ^IJjmc I ' 

yij- 



(65) 
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where stands for 

K = -I'gr"' + - gf"" {{du^;) - (^-^;) O ■ (66) 

As in the case of the toy model, we see that the soft breaking of the BRST symmetry turns out 
to be related to contact terms, enabling us to introduce the extended operator we are looking for. 
Eq.( !65|) expresses precisely the exact invariance of the Gribov-Zwanziger action Sqz under the 
operator 





s-yScz — , 


(67) 


where is given by 


























= 0, 








' yx ' 




= 0, 

— 1 j"-^ 

— ^tM ) 

= . 


(68) 


Like the operator of eq. fl3T| 


, the modified operator is non-local. 


However, unlike s^, it lacks 



nilpotency. Nevertheless, it turns out to be helpful in order to evaluate the vacuum expectation 
value of BRST exact quantities, like the condensate (O | (^^^(a;)^?^* 



3.4 Evaluation of (O | {Tpf{x)ipf{x) - ujf{x)ujf{x)) \ 



To evaluate the condensate | {(p'^^{x)ip'^^{x) — uj'^^{x)uj1^{xy) | O) we start from the condition 

(0|.,(^;^(x)<^(x))|0> = 0. (69) 
Making use of transformations fl68|) . we get 

(o|.,(cj;^(xX(x))|o> =(o|(^;V?-«)|o>+7V^''^/ A (o|A;:(y)(A<-i)^>™^(x)|( 

(70) 

Therefore 

(0|(^?(xX(x) -aj;:''(a;X(x))|0> = -7^^"^ / A (o Al{y) (A1-^)^7 



Al{y){M-X.^T^^) 0) • (71) 



Equation (ITT!) is a non-perturbative Ward identity relating the dimension two condensate to the 
Green's function appearing in the right-hand side. One should also notice that expression (17T|) 
displays explicit dependence from the Gribov parameter 7, as expected from the non-perturbative 
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nature of (O | {^pf^^, 
To the first order approximation 



)|o>- 



i^-X = (z^ 



5*™ + ... , 



(72) 



yx 



SO that 



(0 1 {rj'vt - 1 0) 



d^p 1 



yx 







-^pi^-y){o\A';^{y)^';{x)\o) 



2 „ rabc 



-7 gf 



d^p 1 



(o|A;(p)y.i'=(-p)|o> 



For the propagator ^0 \A'^^{p)Lp^j^{—p) \ O), we have [TT] 



(o|^::(p)¥'^'^(-p)|o> 



7 



V2^?7' • 



(73) 



(74) 



Therefore 



<o|(^V?-a^''0|o> 



3 

3^2 



2 i\_,2 



d^p 1 



(27r)>2 (p4 + iV74) 



(75) 



647r' 



gVN{N' - 1)Y 



sad 



where use has been made of 

j^abc jdbc 

Expression fl75|) agrees with the result already found in [TT] . 



(76) 



3.5 Remarks on the condition (0 |s^O(a:)| 0) = 

In this section we aim at providing a better understanding of the condition (0 |s^0(x)| 0) = 0. To 
that purpose, we shall study the quantity (O ((y'/"*'^A^(a;)cj|f (a;)) | O). 



From eqs. (l68|) . we have 



1 {gr''^ {Dfd') uj';) 1 0) + (0 \gr''A';if'>;\ o> 



Wg'r'' I d% f 



j4„, fmpc 



-l\pb 
yx 



(77) 



Taking into account that expression ^0 | [gf"''"^ (^D^fd'^ uJ'^^ | O) vanishes, due to the discrete sym- 
metry (uJ^^ — > — Uj^'^ , uj^^ of action one gets 



(0 I i^r^^) I 0) ^ (0 \gr^^\ 0) + ,Vf- J A (O \miy) (M-Xx 



(7J 
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Reminding now the expression of the horizon function h{x) of eq. (H7|) , it follows 

(0 \s, (gr'^A'^.uf;) I 0) = (0 0) + f (0 \h{x)\0) = , (79) 

due to equation (!60l) . 

We see thus that the equation (O \ ((y'/'^^'^A^ZUj^'^) [ O) = follows from the gap equation (160|) 
for the Gribov parameter 7. This provides support to the condition (0|s^B(x)|0) = 0, which 
seems to emerge in a natural way within the Gribov-Zwanziger framework. 



3.6 A general formula 

We can now establish a general procedure to evaluate the BRST exact quantity (0 |s0(x)| 0). Let 
us decompose the operator as 

Sj = s + 6^ , (80) 
where s is the BRST operator of eqs.([3]) and 6^ stands for the operator 





= 0, 




= 0, 


5^^{x) 


- 1' 




= 0, 


5,ujf{x) 


= 1^9 




= 0, 




= 0, 




= . 



, ma 
lyx 



IX pa 



yx 



(81) 

Therefore, from the condition 

(0|s^e(a;)|0) = , (82) 

we get 

(0|se(x)|0) = - (O|(5^0(x)|O) . (83) 

We see thus that the characterization of quantities which are BRST exact relies on the operator 
5^ which takes into account the presence of the Gribov horizon through its explicit dependence 
from the parameter 7. 



3.6.1 From condensates to correlation functions 

The formula fl83|) enables us to extend the previous construction to the case of correlation functions 
which are BRST exact. As an example, let us evaluate the two-point correlation function 

- <(xX(y))> = {s {ujf{x)vf{y))) . (84) 
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From the condition 

{s,{ujfix)^';^iy))) = 0, (85) 

we easily get 



4 The refined Gribov-Zwanziger action 

The so called refined Gribov-Zwanziger action was introduced in [Hill]. It takes into account the 
existence of the dimension two condensate | (^^^v^^* — aJ^^a;^'') | O), which is non-vanishing for 
non-vanishing Gribov parameter 7. This condensate reflects the non-trivial dynamics developed 
by the auxiliary flelds (^^^, y'"'', oJ^'', w'^'') ■ This feature has motivated the introduction of the 
refined Gribov-Zwanziger action, in which the operator (jp ) is taken into account 

from the beginning. This is achieved by adding to the Gribov-Zwanziger action the term 

Stp^ = - j d'x /i^ {rM^ - «) , (87) 

which is left invariant by the BRST transformations ([3]). The parameter /i^ is a mass parameter 
which, as much as the Gribov parameter 7, is fixed by a variational principle, see [H IllJ^ . The 
refined Gribov-Zwanziger action is thus defined by 

'S'rgz = 'S'gz + Si^^ . (88) 

Since the term S"^^ is BRST invariant, the action S'rgz displays the same breaking term of the 
Gribov-Zwanziger action S'gz, ^-e. 

-^'g I d'x (c"^D7 [r'^ {y,'; + 1;)) + r'-Ay;) . (89) 

The refined Gribov-Zwanziger action turns out to be renormalizable to all orders, due to the large 
set of Ward identities which can be established [Hill]. It is worth mentioning that the action (188|) 
gives rise to a decoupling type solution for the ghost and gluon propagator. The ghost propagator 
turns out to behave as in the infrared, while the gluon propagator does not vanish at the 
origin in momentum space, exhibiting positivity violatioij^, namely 

The infrared behavior displayed by expression (l90ll is in good agreement with the most recent 
numerical simulations performed on huge lattices [HI [151 [161 [T71 [H] , cls well as with the analytical 

■^The parameter /i is determined by requiring that the Gribov no-pole condition [3] is fulfiUed in the presence 
of the new operator (|87p . so that the Gribov horizon is not crossed. In other words, the introduction of the term 
(j87p can be done in such a way that one remains within the Gribov region [Jl lllj 

■^The positivity violation of the gluon propagator invalidates the interpretation of the gluons as excitations of 
the physical spectrum of the theory, providing a simple understanding of gluon confinement. 
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results found in [191 [20]. Let us underline that, without the introduction of the operator (IHTIl . the 
Gribov-Zwanziger action would give rise to a different behavior for the propagators, namely to an 
enhanced ghost and to a gluon propagator which would vanish at the origin, see also the recent 
discussion presented in 



Let us proceed by showing how the operator s-y generalizes to the present case. From the re- 
fined Gribov-Zwanziger action, we get 



6S, 



RGZ 



5c^ 



abJj 



SO that 



(91) 



he 



1 ^s^ 



RGZ 



((-MV)"' 



da SS 



RGZ 



da 



(92) 



Thus, equation ( ]89ll can be rewritten as 



.^RGz = / d^x (- {DTK). [M-X: ^ - l'9r''Al(x) ((-M V)"') 



RGZ 



xy 



5uj 



(93) 



with A° given in expression ( |66l) . Eq.( l93l) expresses the exact invariance of the refined Gribov- 
Zwanziger action S'rgz under the nonlocal operator s 



(m) 



S^r^^RGZ = 



(94) 



where the operator s^^'* is now given by 



b%x) + J d'y {DZK)y [M-') 
0, 

r:'{x)+^^g I d'y n'A^{y) ((^V) 
, 



_\\ma 
'yx 



yx 



ijj. 







ah 



(95) 



The same calculations performed in the previous sections with the operator within the Gribov- 
Zwanziger action, can be repeated by employing the operator s"^^ and by taking as starting action 
the refined Gribov-Zwanziger action. 
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5 Conclusion 



• In this work we have addressed the issue of the BRST symmetry in presence of the Gribov 
horizon. Our main results are expressed by equations fl67|) — fl68l) . which show that, due to the 
positivity of the Faddeev-Popov operator Ai"-^ within the Gribov region Q, the soft breaking of 
the BRST symmetry exhibited by the Gribov-Zwanziger action can be converted into an exact 
invariance. The resulting modified operator displays an explicit dependence from the Gribov 
parameter 7. We have shown, through few examples, that the operator might be helpful in 
order to evaluate the vacuum expectation value of quantities which are BRST exact. 

As far as the renormalizability properties of the Gribov-Zwanziger action are concerned, the 
operator has not to be considered as a substitute of the BRST operator s. Albeit broken, we 
remind here that the BRST transformations give rise to a set of softly broken Slavnov- Taylor 
identities which are perfectly suitable for a study of the renormalizability of the theory within a 
local framework O [71 [9], [TOl [11], [12] . Rather, the operator has to be seen as a useful device to 
extract non-perturbative Ward identities in order to characterize the vacuum expectation value 
of quantities which are BRST exact. 

• It is worth spending here a few words about the Gribov-Zwanziger (GZ) action and the so- 
called refined Gribov-Zwanziger action (RGZ) introduced in [11 [11]. Both GZ and RGZ models 
are renormalizable and exhibit the same BRST soft breaking. Moreover, as it is apparent from 
the expression (1871) . the GZ and RGZ models differ by a BRST exact term, notably 

irM' - «) = ' P>t) ■ (96) 

As already mentioned, the introduction of this operator in the RGZ action is supported by the 
non- vanishing value of the condensate (O | (^^''v'^^ — tJ^^u;^'') | O), whose non-perturbative effects 
are taken into account from the beginning, through the introduction of expression f[96l) . We point 
out that the operator f[96l) is the only operator which can be added to the GZ action without 
jeopardizing the renormalizability, while keeping the same BRST breaking as well as the same 
field content. In the absence of the Gribov horizon, and thus in the case of an unbroken BRST 
symmetry, the two actions would be physically equivalent. However, in the present case in which 
the BRST symmetry is softly broken, they give rise to rather different predictions, both in four 
and three dimensions. In fact, while the GZ action predicts an infrared enhanced ghost sector 
and a gluon propagator which vanishes at the origin [21], the output of the RGZ action is a ghost 
propagator which is not enhanced in the infrared and a gluon propagator which does not vanish 
at the origin in momentum space. So far, the behavior predicted by the RGZ model seems to be 
in good agreement with the recent numerical data [T^ [T5 | [T6l [T7t [18] . Interestingly, in two dimen- 
sions, both models coalesce [22 |p|, giving an enhanced ghost and a vanishing gluon propagator, in 
agreement with the lattice data in two dimensions [231 [HI [IB]- At present, we are unable to give 
a simple explanation of the different behavior predicted by the GZ and RGZ actions in three and 
four dimensions. We believe that a better understanding of the meaning of the BRST breaking 
and of its deep connection with the Gribov horizon is needed. 



^This is due to fact that in 2d the operator ([M)) cannot be safely introduced [H], due to the appearance of 
infrared divergencies in the Feynman integrals, typical of 2d theories. We notice in fact that in 2d the auxiliary 
fields (^°^, tU"'', w"'') are dimensionless. 
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• Another issue which deserves further investigation is that of the role played by the auxiliary fields 
(jp'^ , (p'^ ,ZJl^ , u'i^'') within the Gribov-Zwanziger framework. As already noticed, the introduction 
of these fields stems from the necessity of localizing the horizon function, eq.([T]), so that a local 
and renormalizable action is obtained from the restriction to the Gribov region Q. These fields 
are thus expected to carry nontrivial information about the Gribov horizon, a fact which can be 
already inferred from the explicit dependence of the dimension two condensate (^^^v'^^ — uj^^ujI^) 
from the Gribov parameter 7, eq.(jl]). Also, as the restriction to the Gribov region deeply modifies 
the infrared behavior of the theory, we believe that a better understanding of the dynamics of 
the auxiliary fields could provide useful information about the confining character of the theory 
as well as about the construction of its physical spectrum. It is worth therefore to add here a 
few remarks on the on-going investigation [21] about the role which these fields might play. For 
illustration purposes, let us rely on an explicit example taken from the Gribov-Zwanziger theory, 
see [6]. Let us start by giving a look at the gluon propagator of the Gribov-Zwanziger theory, 
which can be obtained from the expression ( !90l) by setting = 0, namely 

{Al{m{-k)) = 5- (s,. - ^) ^ , f = , (97) 

As it is apparent, expression fl97|) displays complex poles at k"^ = ±«7^, which invalidates the 
interpretation of gluons as excitations of the physical spectrum. In other words, in the infrared, 
gluons cannot be considered as part of the physical spectrum: they are confined. The physical 
excitations of the theory would correspond to colorless bound states of gluons, i.e. glueballs. We 
expect thus that information about the physical spectrum could be obtained by looking at the 
correlation function of colorless gauge invariant operators as, for example: 

G(rt^/Ae-'M^=WnO)>. (98, 

where F'^{x) = This correlation function is useful in order to investigate the 

properties of the scalar spin zero glueball state. Within the Gribov-Zwanziger theory, it has to 
be evaluated with the Feynman rules obtained by employing the confining gluon propagator (1971) . 
The explicit first order evaluation of expression ( !98ll can be found in [6j, and can be summarized 
as follows: 

G{p) = G^^^'ip) + ^""^^^^(p) . (99) 

The unphysical part, G^'^^^^^{p), displays cuts beginning at the unphysical values = ±4^7^, 
whereas the physical part, G^^^^{p), has a cut beginning at the physical threshold p"^ = —27^. 
Moreover, the spectral function of G^^^^{p) turns out to be positive, so that it obeys the Kallen- 
Lehmann representation [H]. As such, G''^^^'^{p) is an acceptable correlation function for physical 
glueball excitations. What is interesting in expression fl99l) is that a physical cut has emerged in 
the correlation function of gauge invariant quantities, even if it has been evaluated with a gluon 
propagator exhibiting only unphysical complex poles. This is precisely what one expects from a 
confining theory. Gluons are not physical excitations, but the glueball correlation function dis- 
plays a physical singularity. 

Of course, we have still to face the very difficult presence of the unphysical part G'^^^^^^lp). 
Here, we can argue that the presence of the auxiliary fields (^^^, ip'^,uj1^, uj^^ might be a welcome 
feature. In fact, we can figure out that, in the presence of the Gribov horizon, the correlation 
functions of the physical operators should receive contributions from the horizon, which could be 
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accounted for by the auxiliary fields. This reasoning can be made on a quantitative basis by re- 
calling that in ordinary Yang-Mills theories the construction of the physical operators is obtained 
through the cohomology of the nilpotent BRST operator [21 [25]. It is a remarkable fact that, in 
the present case, the BRST operator s of eqs.Q remains nilpotent, so that physical operators 
can be still associated to its cohomology classes. Furthermore, one should observe that the auxil- 
iary fields (ip^^ , (pi;!' ,uJ^ , Lul^^) are assembled in BRST doublets [25], so that they can only appear 
through terms which are BRST exact. Therefore, for the glueball operator Ogiuebix) we would 
write 

Oglueb{x) = F\x) + slZ^{x) , (100) 

where IZ^ is the exact BRST piece which depends on the auxiliary field^. Expression fllOOj) is 
easily seen to be BRST invariant, namely 

S0giueb{x) = . (101) 

For the glueball correlation function, we would get 

{Ogiuebix)Ogi^M) = {F\x)F\y)) + {sk{x,y)) (102) 

where A(x, y) stands for 

A(x,y) = n^{x)F\y) + F\x)n^{y) + n^{x)s{n^{y)) . (103) 

In the absence of the Gribov horizon, i.e. when 7 = 0, the inclusion of the BRST exact term 
sTZry{x) is completely irrelevant, due to the fact the BRST symmetry is unbroken. As a conse- 
quence, {sA{x,y)) = 0, so that the correlator {O giueb{.x)0 giuebiy)) reduces to {F'^{x)F'^{y). How- 
ever, in the presence of the Gribov horizon, i.e. for 7 7^ 0, BRST symmetry is softly broken, 
so that the quantity {sK{x,y)) is now non-vanishing. We could thus search for a suitable term 
TZ^{x) which would enable us to cancel, order by order, the unphysical contributions contained in 
Qnnphjs^^y Thus, the correlation function {Ogiueb{x)Ogiueb{y)) would display only physical poles, 
corresponding to colorless gluon bound states. At the present stage, such a possibility has to be 
regarded as an interesting plausible working hypothesis |23]. Certainly, it is a point worth to be 
investigated. If true, this would mean that the Gribov horizon combines in a nice and friendly 
fashion with the soft breaking of the BRST operator in such a way that the unphyiscal cuts 
present in the correlation functions of physical operators could be removed order by order. The 
presence of the auxiliary fields would be thus a quite welcome feature, as they would give us a 
way to construct the exact BRST term 7ly{x) needed to account for the unphysical cuts. At the 
same time, the existence of the operator would enable us to evaluate the quantity {sA{x, y)) in 
a simple and efficient way. 

• A topic which we also intend to analyse is that of the possible existence of a relationship 
between the non-local operator and the color invariance of the Gribov-Zwanziger action. This 
could open a small window towards a better understanding of the color symmetry in the presence 
of the Gribov horizon. 

• Finally, it is worth mentioning that the construction of the non-local operator outlined 
here can be generalized to the maximal Abelian gauge [271 [2H1 [291 EOl [SH [32] , for which a study 
of the properties of the corresponding Gribov region is available [33] . 

^See also the recent work (261 . 
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